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ABSTRACT 

A  result  of  Strang  and  Fix  states  that  if  the  order  of  controlled 
approximation  from  a  collection  of  locally  supported  elements  is  k,  then 
there  is  a  finite  linear  combination  ft  of  those  elements  and  their 
translates  such  that  any  polynomial  of  degree  less  than  k  can  be  reproduced 


by  Q  and  its  translates.  This  paper  gives  a  counterexample  to  their  result. 
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SIGNIFICANCE  AND  EXPLANATION 


In  finite  element  analysis  it  is  important  to  choose  appropriate  trial 

functions.  One  approach  to  such  a  choice  has  been  described  by  Strang  as 

(  '  '  -11  •'  f  ‘  .  *  r  ‘ 

follows:  Choose  one  or  more  locally  supported  functions  . ..,  rescale 

the  independent  variable,  and  translate  the  functions  just  constructed.  A 
principal  question  in  finite  element  analysis  is  how  well  a  given  function  can 
then  be  approximated  by  those  trial  functions.  This  is  measured  by  the  order 
of  approximation.  Sometimes  it  is  required  that  the  coefficients  in  the 
process  of  approximation  be  bounded.  Such  an  approximation  is  said  to  be 
controlled.  A  result  of  Strang  and  Fix  states  that  if  a  collection  of  locally 
supported  elements  has  controlled  approximation  of  order  k,  then  there  is  a 

>  V. 

finite  linear  combination  ^  of  those  elements  and  their  translates  such  that 

>o  *  ’ 

any  polynomial  of  degree  less  than  k  can  be  reproduced  by  9  and  *ts 

translates.  This  result  has  been  questioned  for  some  time.  Recently,  the 

rapid  development  of  the  theory  of  multivariate  spline  functions  has  given  a 

new  impetus  to  the  controlled  approximation  problem.  It  becomes  significant 

to  answer  the  question  whether  their  result  is  correct.  This  paper  gives  a 

n  y -  -7 :  ,  '<  -  * 

counterexample  to  the  ^result"  of  Strang  and  Fix.  Our-  conclusion  suggests 

- . ...  .;/,_>*  14.  be.  ■ 

that  we  should  redefine  the  controlled  approximation  order^  Hopefully,  the 

result  in  this  paper  will  lead  to  further  development  of  finite  element 

analysis,  j 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  COUNTEREXAMPLE  TO  A  RESULT  JT  STRANG  AND  FIX 
CONCERNING  CONTROLLED  APPROXIMATION 


Ron^Oing  Jia 

This  paper  concerns  the  so-called  controlled  approximation.  This  concept 
was  introduced  by  £S]  in  1970.  In  that  paper,  Strang  described  a  systematic 
approach  to  the  choice  of  trial  functions  used  in  the  finite  element  analysis 
(Also  see  [FS]  for  some  related  results).  To  describe  his  approach,  we  first 
introduce  some  notation.  As  usual,  we  mean  by  W^'Pd^)  the  usual  Sobolev 


space  with  the  norm 


l  lDaul 
|a|<k  Lp 


(Here  we  use  the  standard  multi-index  notation).  The  seminorm 
wk'Pftf1)  is  defined  to  be 


l  »D°ul 


when  k  -  0,  I*  .  is  a  norm,  and  we  write 
k»P 


By  W*'p(aP)  we  denote  the  space  of  all  functions  in  that  have 

c 

compact  support.  By  4^(2°)  we  mean  the  space  of  the  mappings  w  from  Z° 


to  R  which  satisfy 


»wlp  Cl j I w( j) |p)1/p  <  «  . 


Strang's  approach  can  be  described  as  follows:  Choose  one  or  more  functions 
♦)#••• •♦n  in  w£,2(Rn)#  rescale  the  independent  variable,  obtaining 
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h  $^(*/h)  with  h>0  (i  =  1,...,N),  and  translate  the  functions  just 
constructed,  replacing  x  by  x  -  jh.  Thus  we  have  a  family  of  trial 
functions : 

-n/2 

$  .(x)  :=  h  $  -  j),  i  =  1,...,N,  j  e  . 

1,3  in 

When  h  =  1,  we  omit  the  superscript  h  in  .  •  Note  that  the 

i  *  3 

coefficient  h  n/> 2  normalizes  .  in  the  sense  that 

1*3 

1 4>h  I  =  I*  I  . 

*i,j  2  ?i  2 

A  principal  question  in  finite  element  analysis  is  the  degree  of 

approximation  which  can  be  achieved  by  the  span  of  . •  Concerning  this 

i  *3 

problem,  Strang  and  Fix  state  the  following  result  (see  [SF;  Theorem  II]): 

Theorem  A.  Suppose  are  in  W^'2(l^ri).  Then  the  following 

conditions  are  equivalent: 

(i)k  for  each  u  in  wk+1'2  and  h  >  0,  there  are  weights  ,  sue 
*  1*3 


1  •  J 


const  h^+^“8|u|  for  a  <  k  , 

8  JC+ 1  i  4 


Iw" I-  <  const  lul- 
12  2 


(Here  w^1  is  the  mapping  j  +  w*1  j  e  Z11.) 

i  i*3 

(ii)k  there  is  a  finite  linear  combination  ft  of  and  their 

translates  .  such  that 

i*3 

for  |o|  <  k,  l  jaft(x  -  j)  is  a  polynomial  in  x^,...,xn 

jez11 


with  leading  term  x  . 


but  not  ( i ) 


then  we 


If  a  collection  $  =  satisfies  (D^.i  but  not  (i)^* 

say  that  <t  has  controlled  approximation  order  k.  Note  that  (2) 

constrains  the  coefficients  v}  . .  This  is  the  reason  why  [S]  calls  such 

i  >  j 

approximation  controlled. 

The  implication  of  (1)^  from  (ii)k  is  in  no  doubt.  Wiis  has  been 
proved  by  several  authors,  in  particular,  by  [SF]  (see  the  references  cited  by 
[S]  and  [FS2] ) •  But,  the  other  direction,  i.e.,  the  implication  of  (ii)^ 
from  (i)k  has  been  proved  only  for  N  ■  1  by  [SF] »  the  proof  for  the  case 
N  >  1  is  not  discussed  in  [SF] . 

Recently,  the  rapid  development  of  the  theory  of  multivariate  spline 
functions  has  given  a  new  impetus  to  the  controlled  approximation  problem.  In 
[DM] ,  Theorem  A  was  cited  to  determine  the  controlled  approximation  order  from 
certain  bivariate  box  splines.  In  order  to  advance  the  theory  of  controlled 
approximation,  it  becomes  important  to  answer  the  question  whether  (i)^ 
implies  (ii)k* 

As  it  turns  out,  (i)k  does  not  imply  (ii)*  in  general.  We  shall  give 
a  counterexample  in  the  following.  In  our  example,  n  *  2,  N  *4,  and 

«/$4  are  splines  with  compact  support,  which  are  constructed  from  box 
splines.  Thus  we  have  to  recall  box  splines.  Before  doing  so,  we  need  to 
introduce  more  notation.  The  vectors  e^(i“  1,2,3)  in  R^  are  given  by 
ef  *  (1,0),  e2  ■  (0,1)  and  e3  -  (1,1)  • 

The  difference  operators  are  given  by  the  rule 

V  f  ■f-f(.-ei),  for  f  s  Z2  ♦  R  . 

The  differential  operators  are  given  by 


Di  «  <£ —  (i  -  1,2)  and  D3  -  D.,  +  D2  . 


-3- 


Moreover,  we  denote  by  the  space  of  all  polynomials  in  two  variables  of 

total  degree  no  more  than  k.  Let  us  recall  the  definition  of  box  splines 
from  [BH] .  For  a  sequence  =  =  vectors  in  rP,  the  box  spline  M-  is 

defined  to  be  the  distribution  given  by 

n 

<M-,4>  =  /  _♦(  l  AUK JdA 

[0,1] n  i=1 

00 

for  any  C  -locally  supported  function  $.  If  H  consists  of  e1f  e2  and 
e3,  repeated  r,  s  and  t  times,  respectively,  then  we  write  MrfSft 
instead  of  M=.  In  particular, 

Mr,s,0(x1'x2)  *  Mr(x1)Ms(x2)  . 

Here  Mr  is  the  usual  (univariate)  B-spline  of  order  r  at  the  knot  sequence 
0^1^ •• • / r  z 

Mr(x)  -  [o,1,...,r](*  -  x)*  1  , 

where  [pQ ,p1 , . . .p^ ]f  means  the  divided  difference  of  f  at  points 
Pg,P^»  • • • »P  • 

Now  we  can  formulate  our  collection  4  **  {4^  '  0111  can<Jidates 

are  M2>0 M210,  M0f2^  and  M12q.  But  we  want  to  modify  them  to 

suit  our  purpose.  First,  for  the  technical  reason  which  will  be  clear  later, 
we  want  to  shift  those  splines.  Thus  we  define  the  shift  operator  T  by 

Tf  :  *  f(»  +  e3)  . 

Second,  more  important,  those  splines  are  not  in  the  space  W  '  .  So  we  have 
to  smooth  them.  This  can  be  done  by  convolving  them  with  a  sufficiently 
smooth  function.  We  should  be  careful  in  doing  so,  however,  because  we  want 
to  keep  4  from  satisfying  (ii)^.  Thus  we  introduce  the  convolution 
operator  0  as  follows: 

Of  :  *  /f(x1  -  t,x2  -  t)2M1(2t  +  -j)dt 


rCT.VAWW.K,3P^^  .J-  ■.’.>■.«.■  .,■ 

,  **« 

NV> 

> V 

Note  that  0  and 


rtS 

& 

s 

I 

i 


>.-■ 

L\- 
l'  % 

* 


t.y 


m 


. 

c.-. 

p. 

k. 

T- 

£ 

lV 

fm  0 

v.- 

£ 

K 


commute  with  each 

other . 

Let 

0T  ("2,0,1 

)  o 

*2  :» 

0T 

(**2,1,0 

0T  K,2,V 

)  * 

*4  :« 

OT 

O 

«k 

fM 

sT 

Obviously,  ♦,(*  ■  1,...,4)  are  in  W1*2.  Our  result  can  be  stated  in  the 
following 

Theorem.  The  collection  ♦  satisfies  (i)^  but  not  ( 1 i }  ^ . 

Proof.  To  prove  the  first  assertion,  we  start  with  the  box  spline 

This  spline  is  well-known  to  finite  element  analysts  as  the  standard 
linear  element,  which  was  first  introduced  by  Courant  [C] .  The  spline 
^  y  has  the  hexagon 

lX1«1  +  X2®2  +  X3e3'  0  <  Xi  <  ^ 
as  its  support ,  and  ■  1 «  Let 


He  claim 
(3) 


♦  '  "‘Vi.!1  • 


I2p(j)+(*-j)*P  for  any  p  6  v  . 

jea 


To  this  end,  we  recall  the  following  elementary  properties  of  box  splines  (see 
[BH]  ) : 

I  M;C  •  j)  »  1  . 

and 

(4)  oJZaCjjM-t*  -  j))  -  E(Via)(j)M^c  (•  -  j)  if  e±  e  S  . 


For  p  e  *1#  D^p  -  V^p(i  *  1,2)  are  constants,  hence 

Di(ipcj>«1i,y  -  3  *  «3))  -  i<V>(i|M(.1,.2,.3i\.1'-  -  *  +  e3’>  ■  V  • 


-5- 


Moreover , 

£p( j>Mi< 1,1 (0  "  j  +  e3)  *  p(0)  . 

This  shows  that 

Zp(j)(TM,  -  j)  *  p  for  pew,  . 


Let  the  convolution  operator  0  act  on  both  sides  of  this  equality.  We  have 

Zp{ j)(OTM,,  ,)(•  -  j)  ■  op  . 

It  is  easily  seen  that  op  -  p  for  any  pew,.  This  proves  our  claim  (3). 
From  the  foregoing  argument  we  also  have 

Z < D± 4» )  { -  -  j)  -  1  . 

Thus,  by  a  standard  argument  (see  the  proof  of  [BH;  Theorem  6]),  we  have,  for 
u  e  W2'2(«2)D  C2(B2)  and  h  >  0, 


(5)  |u  -  £u(  jh)$(  j-  -  j)|s  2  <  const  h2_S  |  ul  2, 2 ' 

For  simplicity,  we  write 


s-0,1  . 


w(j)  -  hut  jh) 
v  «  Ew(j)^ 

with 


♦5  *  h’Vii  - 1)  • 

o  2 

Since  u  e  L  ,  w  el  for  sufficiently  small  h  >  0,  and 

Iwlg  <  const  lulj  • 

We  want  to  express  v  as  an  infinite  linear  combination  of  < 

2 

is  possible  because  of  (4).  Suppose  y  :  Z  ♦  R  satisfies  V,y  = 
follows  from  (4)  that 


h 


w. 


.  This 
Then  it 


l  w<j>M, ,,,,<*  -  -  l  {V><j>M1,1,1(*  '  j) 

-  0,(1  Y(j)M21  ,(•  -j))  =  (d3  -  d2)Q  y( j)M,(»  -  j)) 

-  i  -  3>  -  l  -  31  * 


-6- 


Recall  that 


♦  ‘  0T(M,,1,1>'  *1  "  CTT<M2.0,1>  and  *2  '  O,<"2,1,0>- 

We  obtain 

l  -  l  <V)(j)*2,j  "  l  (V)(j)*1,j  * 

Thus  v  can  be  expressed  as  an  infinite  linear  combination  of  4  and 

1 1  j 

$2  j •  But  we  have  trouble  in  keeping  the  coefficient  sequence 
2 

V^Y  bounded  in  L  .  To  settle  this  trouble  we  shall  use  Fourier  series* 

From  now  on,  the  letter  i  will  be  reserved  for  the  imaginary  unit  f~-\ . 

2  2  * 

any  a  e  l  (I>  ) ,  we  denote  by  a  the  Fourier  series  given  by 

a(£)  s  -  I  atjje2*1***  ,  5  e  R2, 

where  +  j2?2  i®  the  inner  product  of  j  and  £.  We  have  the 

following  elementary  property: 


(6)  (Vka)  "  S  (aO)  *  a( j  -  ek))e2iri^’^  *  (1  -  exp{2riCk)a  , 


k  • 


Let 


g^S)  :  *  1  -  exp(2*i5k),  k  -  1,2  . 

Then  (6)  tells  us  that 

(7)  (V2a)  »  g2a  -  (g2^g1)(71a)  . 

Let  0  be  the  square  [0,1]  x  [0,1].  By  Parseval's  identity. 


,w,2 ,Q  !  “  (/Qlw|2)  '*m  ,w,2  * 


Thus  if  =  w,  then 


'V'2  *  ,(72Y)  *2,0  "  l<92/gi)  (VY)  *2,Q  *  l(g2/gl,W,2,Q  * 


Since  g2/gi  is  unbounded,  the  norm  Mgj/g^)  w*  2  Q  may  ^  unbounded  in 
general.  Our  scheme  is  to  decompose  w  into  a  sum  of  two  sequences  a  and 


and  for  b 


b.  For  a,  we  seek  a  such  that  V^a  =  a  and  *^2ol2  <  lwl2> 
we  find  0  such  that  =  b  and  <  lwl2«  This  decomposition  of 

can  be  done  in  the  following  way  :  Introduce  the  set 

E  :=  {(C1fC2)  e  q,  |gi(5)|  >  I g2 ( C) | }  , 

and  its  characteristic  function 


y  (E)  ■ =  1  if  £  6  E 

V*'  ’  o  if  5  e  * 


Decompose  w 

into 

wx  + 

£ 

i 

X 

n 

• 

Obviously, 

,;v 

2,0  < 

A 

lwl_  .  and 
2,0 

'“( 1  "  *e»2.Q  <  ,;i2,Q  * 

Hence  we  can 

expand 

A 

«XE 

A 

and  w(1  - 

XE>  in  Fourier  series  : 

wxe<5)  =  £  a(j)e2irij*5  , 


w(t  -  XE>(§)  -  £  b(j)e2,i3’5 


with 


lal„ , Ibl  <  Iwl  ^  =  Iwl, 

2  2  2,0  2 


It  follows  that 

£  w(j)e2lIi^  ^  =  w($)  =  E  (a(j)  +  b(  j  ))e2ni^*^  . 

By  the  uniqueness  of  Fourier  series  expansion,  we  have 

w  =  a  +  b  . 

A  A 

Since  a  =  wxE  ,  and  since  |g2(£)|  <  | g^ ( C ) |  for  £  e  E,  we  have 

A  A 

l<g2/9i>a(Ul  *  la<5)|  for  5  e  q  , 

2  2 

and  therefore  we  can  find  c  e  A  (Z  )  such  that 

A  A 

c  =  (g2/g1)a  . 


-8- 


•  v  •>»  •  •>  a  •  ^  .*»  ^  ,K  n"  ,N  ,>  ^  a  *  . 


-  a  ~  a  *  •  *  •  -  »  a  , 
•»  *i  V.\*  *  *  '• 

*  ‘  *  Os'4  V  *  ■ 

*  a  *  »  .  ■  *  ■  r  V  *•  . 


W 


(8) 


jr-j. 


In  addition. 


.c'2  -  «c«2  Q  <  «a«2  p  »  »a«2  <  iw<2  . 


Furthermore,  (8)  is  equivalent  to 


gtc  *  g2a  ,  or  (7^)  -  (72a)  i 


hence 


7^c  *  72a  . 

Having  obtained  a  and  c,  we  can  find  (uniquely)  a  sequence 
2 

a  :  Z  *  R  such  that  7^a  *  a  and 


]  f  (0,i)  for  j  >  0 
£-1 


a(0,j2) 


for  j2  »  0 


-  I  f(0,£)  for  j,  <  0 


*-j2+1 


71<72a)  «  72 (V^a)  -  72a  -  7^ 


But  (9)  tells  us  that 


72o(0, j2>  *  c(0,j2)  for  all  j2  e  2  I 


therefore 


72o  =  c  . 


Furthermore 


hence 


73<x  *  a  -  a(*  -  e3)  *  a  -  a(»  -  e2>  +  a(»  -  e2)  -  a(»  -  e3) 


72a  +  Vla(*  "  e2}  “  72a  +  a(' 


I73al2  <  2  lal 2  <  2  lw»2 


-  e2 '  1 


t;W, 


»  .  . '‘'.'•".'-  A'.--' • 


^  J 


To  sum  up,  we  can  find  o  :  Z  ♦  R  such  that 


and 


IV2ol2  <  lwl2  ,  »V3al2 
E  a(  j  =  £  (V3a)(  j)^  -  E 


<  2  Iwl 2 


2 

Similarly,  we  can  find  $  :  Z  ♦  R  such  that 

IV1012  <  lwl2  ,  »V30l2  <  2  lwl2 

and 

E  b(j)^  =  E  (736)(j)^j  -  E  (V^H j)f3  . 

Recall  that 

v  «  E  w(j)^  -  I  a(j)^  +  I  b(j)*|j  . 

From  (5)  we  observe  that 

|u  -  vlSf2  <  const  h2“8!u|2  2  ,  s  =  0,1  . 


Thus  we  have  proved  that  the  collection  *  **  ( ^1 '^2 '^3 '^4^  satisfies  (l)j« 
It  remains  to  show  that  *  does  not  satisfy  (ii)^.  This  will  be  done 
by  contradiction.  Suppose  to  the  contrary  that  there  is  a  finite  linear 
combination  ft  of  and  their  translates  ^  such  that 

for  |o|  <  1,  l  j°fl(x  -  j)  is  a  polynomial  in  x^,X2  with 

jez 

leading  term  x°. 


Assume 


‘  4 


(10)  8  -  l  l  c,  _♦«(•  -  m)  , 

A=1  me I  ' 

2 

where  I  is  a  finite  subset  of  Z •  We  want  to  draw  a  contradiction  from  the 
above  assumption.  Let  us  first  take  a  closer  look  to  the  function  By 

straightforward  computation,  we  have 

^.(x.^,)  -  M, ( x.  -  x_  +1)i|»(x.) 


with 


2<V 

5/4) 

for 

-5/4  <  x2 

< 

-  3/4 

i|>(x  )  ■  1 

for 

-3/4  <  ^ 

< 

-1/4 

-2(x2- 

1/4) 

for 

-1/4  <  x2 

< 

-1/4 

0 

for 

x  <  -4/5 

or 

x  > 

This  motivates  us  to  introduce  the  following  functional  F: 

Fg  :  =  lim  [n  g(j  +  e,  j  +  e)  -  D  g(j  ”  E*  T  "  e)l  9  «  W1f"( 

Just  by  computation,  we  obtain 

2  if  j  -  (0,0) 

F^C  -  j)}  -  -2  if  j  -  (1,1) 

0  otherwise  • 

It  follows  that 

f(Z  j2*1C  "  j>)  *  "2  . 

By  (4)  we  have 

1  h  M2,1,0<-  -  J>  -  D3E  (V*2  +  1,/J)M2,1.1(-  -J> 

-  <»,  +  V  £  (j2U2  *  1>''2)»2,)(,l-  -  )>  •  1  22M2,0.11-  '  5)  ' 

hence 

HD  r(l  j2*2C  -  j))  -  f(e  -  j))  -  -2  . 

By  an  argument  similar  to  that  used  in  proving  (5),  we  observe  that 
l  j2$j^*  “  j)  is  in  n1  for  i  -  3  or  4  . 

Therefore 

(12)  F(j>  j2#tC  -  j))  -  0,  4-3  or  4  . 


Now  we  consider  the  function 


On  the  one  hand,  q  e  ir^;  hence  F(q)  =0.  On  the  other  hand,  by  (10), 


l  j,ae  -  j)  -  l  l  <=£,mW  -  m  ' 


m  4=1 


hence  by  (11)  and  (12), 


F(q)  =  l  I  C*"m  W*  *  ®  -5)  *  -2  I  l  c4,m* 
m  4=1  '  4=1  m 


This  shows  that 


l  l  c4,m  =  0 
4=1  m 


Similarly,  one  can  show 


X  l  c4,m  "  0  * 
4=3  m 


Finally,  since  !<£»(•  -  j)  -  1,  we  get 


l  ae  -  j)  -  I  l  ct  (I  ♦*<*  -  ■  -  J))  -  I  l  «  o 

4=1  m  '  j  4=1  m 

This  shows  that  *  does  not  satisfy  (i.i ) <y  •  Our  proof  is  complete. 


«  .**■ •  ««A  *  • 

*>  *••*»»  * •*»* •  ••  •* » • .  •  i 

*  «  •  O  A.'  *  '  %.*  *  '  .  * 


•j'.'.  *« 
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